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After a short review of inlationary preheating, we discuss the development of equi- 
librium in the frameworks of massless A$ 4 model. It is shown that the process is 
characterised by the appearance of Kolmogorov spectra and the evolution towards 
thermal equilibrium follows self-similar dynamics. Simplified kinetic theory gives 
values for all characteristic exponents which are close to what is observed in lattice 
simulations. This allows estimation of the resulting reheating temperature. 

1. Introduction 

The dynamics of equilibration and thermalization of field theories is of 
interest for various reasons. In high-energy physics understanding of these 
processes is crucial for applications to heavy ion collisions and to reheating 
of the early universe after inflation (for review of inflation see 1 ). Inflationary 
problems for which this understanding can be crutial include 

• Baryogenesis. Indeed, generation of baryon asymmetry is possible 
only in a state which is out of equilibrium. Universe is in highly 
non-equilibrium state after preheating. This opens up new possi- 
bilities and mechanisms for baryogenesis 2 . Clearly, the final unswer 
should depend on details of the equilibration process in many "new" 
and "old" (e.g. Affleck-Dine type 3 ) models. 

• The problem of over-abundant gravitino production in supergravity 
models 4 is avoided if reheating temperature is sufficiently low. 

• The abundance of superheavy dark matter and other relics depends 
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on the time of transition from matter to radiation dominated ex- 
pansion during reheating of the universe . 

Thermalization of field theories was discussed in Refs. 6 . However, at present 
the process of thermalization after preheating is still far away from being 
well understood and developed. The problem is that at the preheating stage 
the occupation numbers are very large, of order of the inverse coupling 
constant. In addition, in many models the coherent inflaton oscillations 
does not decay for a long time. Therefore, a simple kinetic approach is not 
applicable. 

Fortunately, the description in terms of classical field theory is valid in 
this situation 7 , and the process of preheating, as well as subsequent ther- 
malization, can be studied on a lattice. In the papers 8 ' 9 this approach 
was adopted. The goal is to integrate the system on a lattice sufficiently 
accurately and sufficiently far in time to be able to see generic features, 
and possibly to the stage, at which the kinetic description becomes a good 
approximation scheme. Several generic rules of thermalization were formu- 
lated in Ref. 8 , like the early equipartition of energy between coupled fields. 
However, the problem is very complicated and there are other unanswered 
important questions, including: what is the final thermalization temper- 
ature, at what stage the kinetic description becomes valid, what is the 
functional form of particle distributions during the thermalization stage, 
etc. These issues were addressed in Ref. 9 in frameworks of the "minimal" 
inflationary model, the massless A3> 4 -theory. 

It was shown that the particle distribution function follows a self-similar 
evolution related to the turbulent transport of wave energy. This allows 
estimation of the reheating temperature, which turns out to be very low 
in the consiered model. The appearence of self-similar regime of evolution 
should be generic for a wide class of models since typical ranges of particle 
momenta at preheating and in thermal equilibrium are widely separated. 

In the first part of this talk a short review of preheating after inflation is 
given. In the second part we report on the recent progress in understanding 
of thermaliztion after preheating. 

2. Preheating 

Inflation solves the flatness and horizon problems of the standard big bang 
cosmology and provides a calculable mechanism for the generation of ini- 
tial density perturbations 1 . At the end of inflation the Universe was in a 
vacuum-like state. In the process of decay of this state and subsequent ther- 
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malization (reheating) all matter content of the universe is created. It was 
realized recently that the initial stage of reheating, dubbed preheating 10 , is 
a very fast process. This initial stage by now is well understood. 10, 11,7,12-17 

2.1. Fast Inflaton Decay 

Bose-stimulation aids the process of creation of bosons. Occupation num- 
bers grow exponentially with time, n = e' 1 ', which results in a fast, ex- 
plosive decay of the inflaton 10 ' 11 and creates large classical fluctuations at 
low momenta for all coupled Bose-fields. This can have a number of ob- 
servable consequences. Examples include non-thermal phase transitions 16 , 
generation of a stochastic background of the gravitational waves 17 , and a 
possibility for a novel mechanism of baryogenesis 2 . The problem is compli- 
cated, but fortunately, the system in this regime of particle creation became 
classical and can be studied on a lattice 7,12,18 . 

On the other hand, occupation numbers of fermions satisfy n < 1 at 
all times because of the Pauli blocking. This may create an impression 
that the fermionic channel of inflaton decay is not important. This is a 
false impression 14 . Production of supreheavy fermions can be more efficient 
compared to bosons 15 . Let us explain why and when this happens. 

2.2. Fermions versus Bosons 

The effective mass of a scalar X (interaction Lagrangian L- mt = \g 2 (j) 2 X 2 ) 
and of a fermion -0 (interaction Lagrangian L mt — 9<H>il>) coupled to the 
inflaton field 4> are given by the following expressions 



The effective mass of a scalar X depends quadratically upon the inflaton 
field strength and therefore effective mass is always larger than the bare 
mass mx ■ In the case of fermions the inflaton field strength enters linearly 
and the effective mass can cross zero. Therefore, superheavy fermions can 
be created during these moments of zero crossing 15 . Indeed, it is easier to 
create a light field and it is the effective mass which counts at creation. 

The Coupling g by itself is not relevant for the process of creation, g 
always comes in combination with the inflaton field strength. To make 
dimcnsionlcss relevant combination we have to rescale g<j> by a typical time 
scale of creation. In the present case this will be the period of inflaton 



scalar X : m 2 ff = m 2 x + g 2 (j) 2 (t) 
fermion ip : m e s = + g<p{t) . 



(1) 
(2) 
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oscillations or inverse inflaton mass. We obtain 

The parameter q determines the strength of particle production caused by 
the oscillations of the inflaton field. It can be very large even when g is 
small since ^ 2 / TO | w 10 12 . 

A fraction of the initial energy density which goes respectively to bosons 
and fermions is plotted in Fig. 1 for several values of q. For fixed q this 
fraction is a function of the ratio of the particle mass, mgn, to the infla- 
ton mass. We see that superheavy fermions are more efficiently created 
compared to bosons, indeed. 




On the other hand, bosons of moderate or small mass (compared to 
inflaton), and sufficiently large coupling (q > 10 4 , which corresponds to 
g 2 > 10 -8 ) are produced explosively. Their density grows exponentially 
with time until back reaction becomes important. Corresponding field vari- 
ances, (X 2 ), can reach large values, larger than in a thermal equilibrium at 
the same energy density leading to interesting physical consequenses which 
were described above. 
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3. Thermalization 

A study of thermalization requires long and accurate integration on a rea- 
sonably large lattice. Therefore we consider the simplest model for reheat- 
ing after inflation: the massless A$ 4 theory. 



3.1. The Model 

With conformal coupling to gravity and after rescaling of the field, ip = $a, 
where a(t) is the cosmological scale factor, the equation of motion in comov- 
ing coordinates describes </? 4 -theory in Minkowski space-time, Dip+Xf 3 = 0. 
At the end of inflation the field is homogeneous, ip = tpo(t). Later on fluc- 
tuations develop, but the homogeneous component of the field which corre- 
sponds to the zero momentum in the Fourier decomposition may be dynami- 
cally important and is referred to as the "zero-mode." In such situations it is 
convenient to make a further rescaling of the field, <f> = (p/ipo(to), and of the 
space-time coordinates, — > v'Av'o^o)^ . This transforms the equation 
of motion into dimensionless and parameter free form, D<f> + <f> 3 = 0. Here 
to corresponds to the initial moment of time, which is the end of inflation 
for us here. In what follows we denote dimensionless time as r. With this 
rescaling the initial condition for the zero-mode oscillations is 4>o(tq) = 1. 
All model dependence on the coupling constant A and on the initial ampli- 
tude of the field oscillations is encoded now in the initial conditions for the 
small (vacuum) fluctuations of the field with non-zero momenta 7 . The phys- 
ical normalization of the inflationary model corresponds to a dimensionful 
initial amplitude of (po(h) ~ 0.3Mpi and a coupling constant A ~ 10~ 13 . 



3.2. Numerical Procedure and Results 

Various quantities were measured on a 3-D cubic lattice and monitored 
both in configuration space (zero mode, <p = (<fi), and the variance, 
var(<j>) = ((j) 2 } — 4>q) and in Fourier space, where the wave amplitudes 
(which correspond to annihilation operators in the quantum problem) were 
defined as 

aW ~ (2^)3/2^ ' W 

The effective frequency uj k = ^Jk 2 + is determined by the effective 

mass n^eff = 3A(0 2 ). Low order correlators containing a^, namely n(k) = 
(a'a), a(k) = (aa) and (aWaa) were measured. The first one, which 
corresponds to the particle occupation numbers, is of prime interest. 
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Occupation numbers at different moments of time are shown in Fig. 2. 
The first stage of the inflaton decay (t < 1500) is characterized by a peaky 
structure of the spectra. The first peak, which corresponds to the paramet- 
ric resonance, is initially at the theoretically predicted value of k <~ 1.27 n . 
Higher peaks appear later (r > 300) because of re-scattering of waves from 
the first resonant peak. With time the first peak moves to the left reflect- 
ing the change in the effective frequency of inflaton oscillations. However, 
after rescaling of particle momenta by the current amplitude of the zero 
mode oscillations, as it is done in Fig. 2, the position of the resonance is 
approximately unchanged. 

The second stage of the evolution (r > 1500) has following characteristic 
features 9 : 

(1) The system is statistically close to Gaussian distribution of field 
amplitudes and conjugated momenta, while phases are random, see 
also 6 ' 8 . 

(2) In accord with this, the "anomalous" correlators are small (but non- 
vanishing), e.g. (Tfe ~ 10~ 2 rife. 

(3) The zero-mode <fio and variance var(4>) are decreasing according to 
the power-laws 4>o ~ t -1 / 3 and var{4>) ~ t~ 2 / 5 . 

(4) The zero-mode is in a non-trivial dynamical equilibrium with the 
bath of highly occupied modes: when the zero-mode is artificially 
removed, it is recreated on a short time-scale (Bose condensation). 

(5) The spectra in the dynamically important region can be described 
by a power low, k~ s with s ~ 1.5—1.7. We see that, despite the fea- 
ture (1) holds, the system is not in a thermal equilibrium. Rather, 
particle distributions correspond to Kolmogorov turbulence. 

(6) The power law is followed by an exponential cut-off at higher k. 
Energy accumulated in particles is concentrated in the region were 
the cut-off starts. It's position is monotously growing, reflecting the 
motion towars thermal equilibrium. 

(7) This motion can be described as a self-similiar evolution 

n(fc,r) = T- q n {kT- p ) , (5) 
the best numerical fits are q w 3.5p and p w 1/5, see Fig. 3. 

3.3. Discussion 

Let us see if these numerical results can be understood analytically. Fea- 
tures (1) and (2) facilitate the use of a simple kinetic approach, rik = Ik, 
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Figure 3. On the right hand side we plot the wave energy per decade found in lattice 
integration at r = 3600, 5100, 7000, 10000. On the left hand side are the same graphs 
transformed according to the relation inverse to Eq. (5). 



where the collision integral for a m-particle interaction is given by 9 

h= fdQ k U k F[n}. 



(6) 
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In d spatial dimensions the integration measure dQk corresponds to m — 1 
integrations over d-dimensional Fourier space. Included in it is the energy- 
momentum conservation <5-functions but not the relativistic l/w(fcj) "on- 
shell" factors, which instead appear in the "matrix element" of the corre- 
sponding process, Uk- This convention makes the discussion of relativistic 
and non-relativistic cases uniform. The function F[n] is a sum of products 
of the type n^. 1 YliLi n ki i where j e {1, . . . , m} with appropriate signs and 
permutations of indices for incoming and outgoing particles. All dynamical 
aspects of turbulence follow from the scaling properties of the system 19,20 . 
Let ujk,nk and Uk have defined weights under a £-rescaling of the Fourier 
space, 



The weight of the full collision integral under this reparametrization is 



Then stationary turbulence with constant energy flux over momentum 
space is characterised by a power-law distribution function, nk ~ k~ s , 
where s — d + (3/(m — 1). The scaling properties also give the exponents 
of the self-similar distribution, Eq. (5). A system with energy conservation 
in particles and £ = t~ p has q = Ap and p = l/((m — l)a — (3), while for 
stationary turbulence exponent p should be (m — 1) times larger. 

For a massless A</> 4 -theory in three spatial dimensions and four-particle 
interaction we have m = 4, /? = —4a and a = 1. In this case s = 5/3 
and p = 1/7 9 . Obtained value for p is somewhat smaller then observed 
numerically, p = 1/5. For three-particle interaction (the fourth particle 
belongs to the condensate in this case and the matrix clement contains an 
additional factor cf) ) the stationary turbulence is characterized by s = 3/2 
and p has even smaller value compared to 4-particle interaction. Numerical 
results 9 do not allow to distinguish between 5/3 and 3/2, s rather fluctuates 
between these two numbers, while 1/7 for p gives a fit to the data not 
as good as displayed in Fig. 3. However, taking into account the non- 
stanionarity of the source (influx of energy from the zero-mode is still non- 
negligible) we obtain 3.5p and p w 1/6 9 . This should be considered as 
satisfactory agreement between numerical experiment and a theory, given 
the simplifications which were made. 



^U{k u ...,k m ) 



(7) 




(8) 
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3.4. Equilibration time and temperature 

At late times, with decreased influence of the zero-mode, our simple an- 
alytical description should become more precise. On the other hand we 
can expect that the self-similarity persists. In this case solution Eq. (5) 
with p = 1/7, which determines the transport of energy to higher mo- 
menta, should be valid. In classical frameworks it will be valid at all later 
times. Therefore, the time needed to reach equilibrium can be defined 
as the time when classical approximation breaks down. In other words, 
equilibrium will be reached after occupation numbers in the region of the 
peak, see Fig. 3, will become of order one. Using Eq. 5 we find 9 for the 
equilibration time r ~ A~ 7//4 ~ 10 23 , where in the second equality the 
normalization to the inflationary model is assumed. Rotating back from 
the conformal reference frame, we obtain for the reheating temperature 
T R ~ k max /a(T) ~ AVo(to) ~ lCr 26 M P1 - 100 eV 9 . This result coin- 
sides with what could have been obtained in "naive" perturbation theory. 
Namely, equating the rate of scattering in thermal equilibrium to the Hub- 
ble expansion rate we get in this model T ~ A 2 Mpi. 

4. Conclusions 

Reheating after preheating appears to be a rather slow process. Although 
the "effective temperature" measured at low momentum modes during pre- 
heating may be high, in the model we have considered the resulting true 
temperature is parametrically the same as what could have been obtained 
in "naive" perturbation theory. We anticipate this result should be ap- 
plicable to more realistic models of inflation. Note that realistic models 
involve many fields and interactions and the largest coupling constants will 
determine the true temperature. 
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